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1. The present paper is suggested by that of Dr. H. F. Baker in the 'Proceedings 
of the London Mathematical Society/ vol. xxxv., p. 333, "On the Integration of 
Linear Differential Equations." In that paper a linear ordinary differential equation 
of order n is considered as derived from a system of n linear simultaneous differential 
equations 

or, in abbreviated notation, 

dxjdt =• ux, 
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where u is a square matrix of n rows and columns whose elements are functions 
of t, and x denotes a column of n independent variables. 

A symbolic solution of this system is there given and denoted by the symbol Cl(u). 
This is a matrix of n rows and columns formed from u as follows : — Q(<f>) is the matrix 
of which each element is the ^-integral from t to t of the corresponding element of <£, 
<j> being any matrix of n rows and columns ; then 

a(u) — l-i-Qu J rQuQ^ + QuQuQu t ..ad inf., 

where the operator Q affects the whole of the part following it in any term. 
Each column of this matrix Cl(u) gives a set of solutions of the equations 

and since 0,(u) = 1 for t = t , these n sets are linearly independent, so that Q,(u) may 
be considered as a complete solution of the system. 

Part II. of the same paper discusses the form of the matrix Sl(u) in the neighbour- 
hood of a point at which the elements of the matrix %i have poles of the first order, 
or in the neighbourhood of which the integrals of the original equation are all 
^ regular T 

It is there shown that if t = be such a point, a matrix 



<t> 



i 
t 



Bu 



A-»» 



c*(t/to) ti - 9 ; c 31 (t/t y 

• • • • Ut 
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can be found, in which all elements to the left of the diagonal are zero, in which 
Cji = unless 6 { — 6j is zero or a positive integer, such that il(u) is of the form 

Gf2(</->)G "~\ 

where G is a matrix whose elements are converging power series in t 9 and G is the 
value of G at t = t . 

The form of <j) is such as to put in evidence what are known as Hamburger's sub- 
groups of integrals associated with the fundamental equation of the singularity ; the 
method is, in fact, a means of analysing the matrix n(u) into a product of matrices, 
of which one is expressible in finite terms and shows the nature of the point as a 
singularity of the solution. 

The object of the following investigation is to see how far, under what conditions, 
and in what form, such an analysis can be effected for equations having poles of a 
higher order than unity in the elements of the matrix u. 

It is known that if in the neighbourhood of infinity the equation is of the form 

1 1 ] 1 z in ~ r) = 0, 



»-i 



2 W + 2 j_p r (aj) + P r (- 
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p r being a polynomial of degree pr 9 and V r (l/x) a series of positive integral powers 
of l/x, the equation has a set of formal solutions of the form 







e Or X <Tr g gjg^ r _ 1 ### ^ ? 



— 00 



where O r is a polynomial of degree p+1, provided a certain determinantal equation 
has its roots all different. 

The case in which these roots are not all different is discussed by Fabry (' These, 
Faculty des Sciences, Paris/ 1885), where he introduces the so-called Subnormal 
Integrals, viz., integrals of the above form in a variable x l!k , h being a positive integer. 

The investigation carried out in the following bears the same relation to the 
discussion of these normal and subnormal integrals that Part II. of the paper quoted 
at the outset bears to the ordinary analysis of the integrals of an equation in the 
neighbourhood of a point near which all the integrals are regular. 

2. Throughout the discussion the neighbourhood of the point t = will be under 
consideration, the coefficient p r being supposed to have a pole of order toy at this 
point. 

Let p + 1 be the least integer not less than the greatest of the quantities m r /r. 
The equation may then be considered as belonging to the more general type 

where P r (0 is holomorphic near t = 0. 

This equation may be reduced to a linear system of simultaneous equations as 
follows (vide < Proc. Lond. Math. Soc./ vol. xxxv., p. 344) : — 

Put x x = z, x 2 = t p+ W\ ... x r+1 = t r{p+1 ¥ r \ r = 1, ..., n-1. 

The n equations then satisfy the system of n equations 



dx 
dt 



r 



j 



0. 



0. 
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o, 



fp + p 



o, 



0. 
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/2? + l> 



(n-2)£±! 



j^p+l 



Sow 
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L+( n -l)£±l 

5 ; 



where Qi*..Q n are series of positive integral powers of t. This system belongs to the 
more general form 
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where ^ +1 .../3 1 ... are square matrices of constants. 

The most general equation of this form will be considered. 

If /a be any matrix of constants and y = fix, the n quantities y satisfy the system 

_ i CO'if (a/'Ju / i \ 

^ dt = ^" ==waj=== ^ )^> 
or 

C £ = ( W" 1 ) y. 

Let /x be now chosen so that (/xop+j/x" 1 ) is of canonical form as follows : — (i.) It has 
zero everywhere save in the diagonal and the n—1 places immediately to the right of 
it; (ii.) The diagonal consists of the roots of the equation \a. p+1 — p\ = 0, equal roots 
occupying consecutive places; (iii.) The elements to the right of the diagonal consist 
of (e 1 —l) unities, then a zero, (e 2 — -1) unities, a zero, and so on (' Proc. Lond. Math. 
Soc./ vol. xxxv., p. 352). 

Form now the matrices (/xa^/x" 1 )...^^.^" 1 ) ; the equation is then replaced by an 
equation of exactly similar form, the matrices a p ... being still any matrices whatever, 
but a p+1 being of the canonical form. 

3. The equation being denoted by 

dyjdt = uy 9 
if i) be any. solution of the equation 

(A) drj/dt = UTj—rix, 

X being an arbitrary matrix, we have 

d 

— { v a( x )} ^vx n (x) + ( u v~vx) n x 

= u{ v n( x )}, 

so that r)Q (x) ' m a matrix satisfying the equation in question, 

hi what follows we are concerned with the form of a solution more than the actual 
convergence and existence of the same. It is therefore important to notice that if t? 
be a diverging power series formally satisfying equation (A), yjfl (x) ma J t> e still 
considered as a formal solution of the original equation, the only condition necessary 
to secure its actual existence being the convergence of 77. 

If 7] be convergent, the solution may be particularized by adding the factor r) ~~ l , 
i.e., rjll (x) rjo" 1 is the solution reducing to unity at t = t Q . 

The main investigation to be carried out is that of a simple form for the matrix ^, 
such that the subsidiary equation (A) may have a formal solution in the form of a 
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matrix whose elements are series of positive integral powers of t, reducing for t = 
to the matrix unity. 

4. Owing to the much greater simplicity of the case in which the equation 
| a P +i—p\ = has all its roots different, it will be treated first separately. The result 
obtained is as follows : — 

A matrix ^ can be determined uniquely of the form 



t p+1 t p t 

where ^ +1 ...^ are matrices of constants in which all elements save those in the 
diagonal are zero, such that there is a formal solution 



V a (x)> 

where the matrix rj is made up of series of positive integral powers of t 
diverging — and reducing for t ■ = to the matrix unity. 
Consider the equation 



generally 



b 



(B) 



drj 



o 



irj—rj 






)• = ! 
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e r n 



The roots of | ot p+l —p 



— being unequal, the matrix a p+1 will have zero elements 
except in the diagonal; the diagonal elements will be p u p 2 , ...p„, the roots of the 
equation. 

If the equation (B) is satisfied by the matrix 

where x 9 y, z . . = denote columns of elements of the form 

££ ^ X'q "T" X-yV ~T" X^ i • • • j 

the coefficients £c r , «/ r ... being columns of constants x r °, x*, x r 2 , &c, these constants 
satisfy the following equations : — 



X. 



\ a p+l C7 p + ij Xq — 0, 

\ a p+l~~~V p+l) X X ~Y \CL p U p } X = u, 
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\ a P+i " p+i)Wp~)~ ••• +( a i — ^i/^o = 0, 



e « 



A precisely similar set of equations gives the relations connecting the constants 
y and 6 2 p+1 ...6 l 2 . 

The equations X just written determine uniquely a set of values for 6 1 p+l >,.0 1 1 and 
the coefficients x . . . . 

The first of these equations gives 

Since x is to be equal to : we must have l p+1 = p l9 and these equations are 

W 
then satisfied. , N 

gives # 2 ^ +1 = p 2 ; and so for 



Similarly the first of the y equations with y = 



1 

0/ 



the other columns. 

The second of the equations X written more fully gives 

a p U —0p 1 = 0> (pr— pi)Wi r + a p lr = 0> r = 2, ..., n. 

These then determine x 1 save for its first element, in place of which a unique value 
is given for Op 1 . 

The third equation X in full gives 

OLp X\ ~r • * • i Ot-p Xi ~j~ Gt p—i v p—\ ~— \) j 

( p r "~ Pi) ^2 • ^» *^1 • ^p *^1 • • • I ( #p t/^ )^l *T" . . • "i O^ ^x T" 0t 2j_i ^ U, r :==z ( Zi) . . . , ?? ). 

Of these, the first gives #^-1, while the following determine x 2 save for its first 
element, but only in terms of the yet undetermined x 1 . 
Of the next group, the first equation is 

Otp X2 • • ♦ 1 &>p X2 -T\GL p—i"-" V p—i)Xi ~\- OL p—iXi -T . • . ~r Ot p_iXi "T OL p—c^—V p~2 == ^« 

This equation apparently involves the unknown x± explicitly, and also, in x 2 2 ...x 2 n , 
implicitly. 

But the whole coefficient of x^ is 



n / Ik \ n 

2 V Pk-pj 2 



n n 

:= 2i0tp X]_ "~~ JLfOLp Jt\ —~- U, 
2 2 



so that 6 l p- 2 is given independently of a^\ 



SATISFYING LINEAR DIFFERENTIAL EQUATIONS. 7 

The remaining equations of this group give x% except for x 3 l , in terms of x x x and x 2 \ 

Proceeding in this way as far as the (jP+l) fch equation of X, x l ...x p are all found 
except for their first elements, while the first elements of these equations give 
8 l P +i--@i l > @ l P +r being not a priori independent of Xi...x* hitherto unknown. 

It has been shown above, however, that the determination of 6 l p ^. 1 does not require 
a knowledge of x x l . 

In general, in fact, 6 1 p ~ r is given independently of x±...x r l . 

To prove this, the way in which x^ enters into x k r+1 will be first considered. This 
may be stated as follows : — 

The coefficient of x x l in x k r+l is equal to that part of x r k which is independent of 

For r = 1 this is at once seen by writing down the equations 

(p a -p 1 )x l k +u p 1 * = 0, 

\Pfc" mm ' Pi) ^2 » * ®"p *^1 ' • • • "» ^p X\ ~T CC p„j ~ (J. 

In general the equation for x k r+l is 



7* 



\pk Pi)'*' r+1 — ^ \ a p-r+sfts "• * • • •" \ a p-r+s " p~r+s)^s ...+#. p - r+s X s j -f- CL p ^ r . 



s=l 



Assuming the statement above to be true for 1, 2...r, and that 0p, ..., 0p- r+1 are 
independent of ^ 1 1 ...^ 1 r -~u the above equation shows that the coefficient of x± in 
— (p k —p 1 )x k r+l is the part independent of asi 1 ,..., < x / * Ln 



.9 = 1 

i.e., in 

so that under the above assumptions the statement holds for 1, 2..., r+1. 
Also, under the same assumptions, from the equation giving 6 l p - r , viz., 

_, 21 sy, 2 1 1 nl^y, w 

-\-\(X p—i' U p—ijX r ~-i ~T" & p -~\X r _i ~T . . * T* OC p—iX r —i "T • • • 

-h^Ot ^„ r + 1 V p-. r + i)Xi ~T . . . ~T Gt ■ ? ;~r + r^'l T" (^ ,p_ r U p — r J = U, 

we deduce that the coefficient of ^x 1 in ^_ r is the part independent of x x l in 

* 

21^2 ■ 1 nl,y,n 

T" 0t p „iX r _2 "T* . . t "T* Ct ^_2^3 r _2 "T" • « • 

This expression differs only from the left-hand member of the equation for l p ~ r+1 
by multiples of x^.^x 1 ^^ and therefore, on the assumption that this equation gives 
d x p - x independently of #/..., the part independent of these quantities in the above 
xpression must, when 0\- r + x is determined, vanish, so that d\- r is independent of x?. 



0k 
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Now the way in which the successive equations follow one another shows that the 
coefficient of x^ in k l is equal to that of x r l in k+r -i. 

Thus $V~ r+ i being independent of x^ 6 1 p ^ r is independent of x 2 l , and in general, 
O l v ^ r+k (lc = l...r— 1) being all independent of x ± l , l p „ T does not contain Xi„.x l r ^ 

Thus, if the assumptions made on p. 7 are satisfied for any particular value of r 
less than p, they are satisfied for a value of r one greater than that value. 

For r = 1 the statements have been justified, and it follows therefore that ^ 1 1?+1 ...^ 1 1 
are all determined uniquely without the knowledge of x^ 9 x 2 1 >.* from the first (jp+1) 
of the equations X, and by the same equations x x .„x p are found, except for their first 
elements, the expressions obtained containing those first elements. 

5. Consider now the (p + 2) th equation X in regard to its first element. 

As before, this will be independent of x p 1 .*.x 2 l ; but on account of the extra term 
arising from dyjfdt, which now enters for the first time, the coefficient of x x l is not 
zero. It is, in fact, —1. 

Thus, the quantities ^ 1 1 ...^ +1 being now known, this equation gives x x \ 

Similarly, the next group's first member will contain the term —2x 2 l but will not 
contain ^ 3 1 ...,x 1 i?+1 , and will therefore give x 2 after x^ is found. 

Thus all the elements x£ are determined successively, and returning to the 
expressions for x k r (r>l) in terms of these and substituting the values so found, all 
these are given also. 

The equations for the columns y 9 %, &c, being treated in the same way, give the 
corresponding ffs uniquely, and also the coefficients in the series of which these 
columns are composed. 

Thus it is shown that when the " characteristic equation" |a p +i— />[ = has its 

roots all different, the equation 

dy/dt = uy y 

where 

a p+l being in its canonical form, possesses a unique formal solution in the form 



y a \J0T + p>~ + *** + £ ) 



where the elements of Xz>--Xi 110 ^ * n ^e diagonal are zero, and the elements of ij are 
power series in t, reducing for t = to the matrix unity. 

The matrix 0(^±r + ... + &) can at once be written in the form a>/w , where w is 

\t p+1 t) 

a matrix whose non-diagonal elements are zero, and whose Jc th diagonal element is 






and oi is the value of a) at t = f . 
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If the series rj happen to be convergent, the solution which reduces to unity at 
t = t can at once be written in the form tjojoiq" 1 ^ 1 . 

Applied to the system formed from a particular linear ordinary equation we have 
at once the result referred to on p. 3 (v. Schlesinger, ' Lin. Diff.-Gleichungen,' 
vol. L, pp. 341 ff. # ). 

As a simple example of the application of the method we ma,y take the well-known equation 



d 2 w _ a -f* ftz + yz 2 
dz 2 # 4 



w. 



Putting W\ — w and w% = z 2 dwjdt 

dwi/dt 



w 2 



which in matrix notation is 



The equation 



-P 1 
a — p 

'a 1 



, dw 2 /dt = ( % + £ + y Wi + - w 2 , 
\z 2 z ) z 

= gives /) = ± Va, and the matrix /x needed to transform the first matrix 



to canonical form is ( ^ * ) , so that the equation is transformed to 



dujdt — 



i/v 



CL 



& 







\ 1 1 2 ^ 







2 J 



a. 



+ 



y 



i 



\/V * \ 1 + -A. i - _JL / 2 V a V 1 

2 V a 2 ^ 



1 



w. 



Vu, say. 



The subsidiary equations are 

drjfdt — Yr] - 7) 



V J* 

zA o 



\ + i/i+i\ 



o \i 



- ^a/ £ \ 0, 1 ~ i^- 



, where X = -£-=, 



which give as the general relations connecting the coefficients of the first column of rj, putting p 
2=1- JX, 



1 + R 



%& 



Hence 

and therefore 









r 



(L). 



- 2 s/ocx 2 n+ i + (p + n) (xn 1 - ^ 2 ) = 0, 
-2 A ya^i 2 + (i? + %--l)(aJ :i w -.i-a; 2 n- i) = ........ 



(**•)> 



which with the first equation gives 



2 Va . n (xn 1 - x n 2 ) = (x\_i - aVi){y + (2 - ft) (# + % -• 1 } . 

Thus a recurrence formula is established for the quantities x^ -x n 2 in terms of which x\ l+ i and x 2 n+ i 
can be at once expressed. 



* With reference to Schlesinger's demonstration of this result, see a note by the author in the 
'Messenger of Mathematics/ January, 1905. 
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The series for x will terminate if for any value of n 

, . c y + (£-*&) (?+**- i) = o> 

i.e., if 

# -7 + (JX + w)(JX + w-l) = 0, 

a.e., if 

(A + 2^~l) 2 -(4y + l) = 
for some value of n. 

The series for y will similarly terminate if 

(-A-f2m~l) 2 -(4y + l) 
vanish for some value of n. 

Both these are certainly satisfied if 

A = q and 4y 4- 1 = ^ 2 , 
where p, ^ are any integers of which one is odd and the other even. 

6. We pass now to the case where the characteristic equation \a p+1 — p\ = has its 
roots not all unequal, and the analysis becomes a good deal more intricate with the 
less simple canonical form of the matrix a p+1 as stated on p. 4. It will be remembered 
that the numbers e l3 e 2 ... there used are the powers of (pi~~p) in the elementary- 
divisors of \ot p+1 ~ p\ with respect to the root p x of this equation of multiplicity L In 
the case of the system obtained on p. 4 from a single equation of order n> we may 
prove that e x = l 9 e 2 = e 3 = . . . = 0, 

For the matrix (a p+1 —p) is of the form 




'0 1 


. 


. . 





1 . 


. . 


« • 

. 


• 
• 


> * • • 

..01 


{a b 


• 


. . . h 



The minor of the quantity " k" in the determinant of this matrix is simply 
(— /o) re "~\ Thus the elementary divisors are certainly merely unity with respect to any 
non-zero roots. 

If there be a multiple zero root, however, since the minor of "a" is unity, the 
elementary divisors with respect to this root are all simply p°. 

Thus for such a system we have for each multiple root e 2 = e 3 . . . = ; so that in the 
canonical form of ot p+1 , if 

n — n /v l ' +1 ' * — 1 

Pi — Pi+li ^p+l — x ? 

and if 

^ -4- ^ „i+l.i A 

Pi -^ Pi+U a p+l — u « 

Such systems being by far the most important in practice and also considerably 
simpler to work out, the full discussion will be restricted to systems of this type. It 
may be pointed out that the most general system can be solved by means of the 
solution of systems of this more restricted type, for from the general system 
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OCX IL i • ••/!• j • (* i it 

= — ^ x, where u is a power series in t, a linear equation ot order n and rank p 

near t = can be obtained for each row of the matrix x, and this equation can be 
solved by the solution of a linear system of the restricted type in question. 

Of the matrix x to be used here, the following properties will be presupposed : — 

(i.) It is to be of the form ^ E ^+^ + ...-i-^- 1 > where each of the matrices 

i? t p t 

Xi-"Xp+i ^ as a ^ elements to the left of the diagonal zero, 
(ii.) The diagonal elements of these matrices are to be numerical constants 

denoted as before by $ r s (r = 1, ... , p + 1 ; s = 1, ... , n). 
(iii.) All the other non-zero elements of XiXz'-Xp+i are ^° be constants, while 

the other elements of ^i m &y contain t, but only to positive integral 

powers (cf. the matrix x m Dr. Baker's paper, loo. cit.), 

8. As before, the matrix rj, which is a formal solution of the subsidiary equation 

drj/dt = urj—rjx, 
will be supposed to be formed of the columns 

X == Xq ~T" Xyj ~r • • • 5 

y = ^o+yi<+—, 

and the equations for the coefficients x r y r ... are the same as the equations X (p. 5). 

But the detailed form of these equations is quite different. The first of them 
(oL p+1 — d l p+1 )x = is still satisfied by 

\!/ 

Supposing now p 3 to be a root of multiplicity e b the second equation X is in full 

X\ ~r ®"p "~" £/« == U , 
^i 3 + V 2 = °> 



• • 



< e) = o, 



♦ •»•••»!» 



( / o n+e2 - / ) 1 )^ + ^+a ?J t ' e >-' e2 = 0, 



where 

C 2 
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These equations manifestly determine x x except for its first and second elements, 
the second being known as soon as p l Is. 

We are also at once faced with a condition necessary for the possibility of the 
solution under the assumptions made as to the form of ^, viz. : — 

This condition arises from the e^ equation of the set, and as, in the ensuing 

discussion, the e/' 31 equation of each set Is most important we shall here introduce a 

notation for it, viz., X r will stand for the e^ 11 equation of the (r+l) th set; i.e., of the 

set 

( a jo+i " p+i) x r -\- •«• == ^* 

This equation will not contain any element of % r . 

A similar notation will be adopted for the equations Y, Z. . • for the coefficients in 
the other columns of rj. 

If the second element of the first row of x ?J +i he c 2l9 the equations Y are 

• *»•#.♦•*«*»» 

Of these the first is satisfied by 




2A> — a 1 3 ^ p+i — P2 ~~ pu 

provided we take e 21 = 1 = corresponding element In ap +1 . 

Considering each of the columns in succession we have thus, with tj q = 1, 

The second of the equations Y gives 

yi 2 +*p 21 -tti 1 = °> 

• •••••*• 



which, when ^ is known, determine ^ save for its first element, and its third until 
P 2 is known. 

The exceptional equation Y u ol^—x^ = 0, gives us again a necessary condition for 
the possibility of the solution in view, o/^ + e^ 1 ' 63 "*"" 1 = 0. 
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Similarly the equation Z x gives 

and so on for the equations for each of the first e x columns of rj. 

For the (e! + l) th column, however, the non-diagonal term of Xp+i = and the 
equations for this column do not contain the elements of the preceding column. 

In fact, the e 2 columns beginning with the (€!+l) th form a group related to one 
another in just the same manner as the first e x are. We obtain from them, as from 
the latter, the conditions 

a/ 1+1 ' ei+e2 "* 1 + a/ 1+2 ' ei+62 = 0, 



a/ li ' 1 ' ei+1 +...+a/ 1+e2,ei+e2 = °> 
and so for the columns associated with each group of equal roots. These and other 
conditions which arise in the course of the investigation will be called " equations of 
condition." Supposing those already found to be satisfied, we may return to the 
solution of the equations X, and of these the following statements are to be 
proved : — 

I. The first Cx—l equations of the (r+l) th set determine x r 2 ...x r €i in terms of 
^ 1 1 ...^ 1 r -i an( i l p+i'"0 l p-r+ii the equations from the (e 1 + l) th onwards give x r * 1+1 . . .x r n 
in terms of the same quantities. 

II. When the values thus found are substituted in the equation X r+1 , the resulting 
left-hand member is independent of the undetermined quantities # 1 1 *»«av L > p 1 ...0 l p - r for 
all values of r up to (e 2 — 2), but for r = e x — 1 it is independent of all save P ; in fact, 
the equation X ei is an algebraic equation of degree e x for 9 P and contains no other 
undetermined quantity. 

III. Supposing one root of this to be chosen for the value of P > and the equations Y 
to be treated in the same way, Y ei _ x will be an equation of degree (e x — 1) in #/, 
whose roots are exactly the remaining roots of X ei . 

IV. Similarly Z e] _ 2 furnishes an equation of degree ex— 2 for #/, whose roots are the 
remaining roots of Y ei _ 1? and therefore of X ei , and so on. 

Thus 6p ...dp 1 are the roots of the equation X ei . 

V. The values of 6 P ^... subsequently obtained in association with each of these 
roots will be the same in whatever order they are taken. 

Of these I. does not require proof. 

With regard to II., the proof that the equations do not contain the undetermined 
x's follows exactly the same lines as the corresponding proof when the characteristic 
equation has its roots all different (vide pp. 7-8). 

The proof that they do not contain (9/... until r = e 1 — 1 requires considerations of a 
different kind involving the equations X, Y, Z... simultaneously. 
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9. Consider the system of equations X 1 derived from X by changing 9 p l into 9 
and x into x 1 , viz. : — 



2 



and let these be treated in exactly the same way as the equations X, the undetermined 
elements of x 1 being supposed the same as those of x. 

From the two sets of equations X, X 1 let a new set be formed by subtracting 
corresponding members of X and X 1 and dividing each remainder by d p —6 2 > and let 
this new system be denoted by 

& p A. = i__/i 2 == U * 

The expressions for a; 1 obtained from X 1 in terms of 0^+1* •• will be identical with 
those obtained from X in terms of the same with 6 p changed into P 2 . 

/y» ... /"v* 



Let A p x denote the expression -^ — ^ • 

p p 

Then A^X is the system of equations 

= 0, (ot^+i u p+ iJA p Xi x = 0, 

( a ^+i~-c/ p+i)ApX 2 + \tt> p u p )A p x l —Xi = 0, 

• • • • • « • » S » • • » J 



Further A p (x x ) = f = 2/o? & n d 0^+! = 6^ +1 




Thus these equations are identical in form with the equations Y, except that 6 2 p - k 
is replaced by d l p ~. k , &>0. Thus if from Y the %/s be calculated as in L, p. 13, and 
from A^X the quantities A p x be similarly determined, the only difference between 
Aa? r +i and y r will be in the substitution of x±... for 2/1 1 ... and 1 p ~ k for 6 2 p - k , &>0, and 
(0 1 1 + l)for^ 1 2 . 

Thus if we substitute the values of A p x thus found in A^X r+1 the result will differ 
from Y r only by the same substitutions. 

In a similar way, denoting by A^Y the difference equations 

Y-Y 1 

Wp Vp 

A^Y will differ from Z r _ x only by the substitution of yi 1 ... for %^ y 9 2 V ~- Y for fl 3 ^..., and 
6 2 + 1 for di> and so for the remainder of the first e x columns. 
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In exactly the same way if X 2 denote the system derived from X by the change of 
6 1 p -i into 6 2 p - u and the equations 

—A 



f) 1 — ft 2 



be formed, the quantities A J ,_ 1 aj r+2 differ from y r only by the substitution of 6 P for p 2 9 
6 l p „ 2 for 6 2 p - 2 ... and 6^ + 2 for ft 2 . The same operator A p _i applied to the equations Y 
connects the columns y and #, and so on. 

Similarly, operators A^_ 2 ...A 2 may be defined. 

X— -X^ 
Lastly, an operator A x may be defined so that the equation A X X is ^ — ^ , where 

X. ip) denotes the equations X with O^—p substituted for X \ Then the equation Y r 
will, when y*,... are replaced by Xi 9 ... and #/,...., 2 2 by #/, ..., 2 x , become the 
equation A^X^ +r . 

Consider now the equation Z x ; it is independent of #/, %\..., and therefore reduces 
to a simple numerical constant which must be zero (p. 13). 

But Y is a polynomial in P . It can therefore be only of the first degree, since 
A^Y 2 is independent of P \ is, in fact, the same as Z l5 viz., zero, so that Y 2 does not 
contain 6 P . It must then, like Z u be only a constant, and must therefore vanish 
identically, so that Y 2 = gives a further " equation of condition." 

Hence again X 3 cannot contain P , and the operator A p - X connecting it with Y x 
shows that it cannot contain 0\-i. Thus X 3 again must be a vanishing constant, 
giving another "equation of condition." 

Similarly, starting with the corresponding equation of the fourth column, we 
find A^X 4 = 0, so that X 4 must be independent of P . Also A^-iX* will vanish 
identically, so that X 4 is independent of 6 l p - u and similarly it is independent of 
ft 1 

v p ~2* • •• 

Thus if €x > 4, X 4 reduces to a mere constant which, as before, must vanish. 

The process may clearly be carried on as far as the €i h column, so that the equations 
Xx.-.X^-x all give equations of condition, as do also Y^^Y^, Z^.Z^, etc. Starting 
now from the second equation of the e x th column, 



th 



where f denotes the €i h column of 77, it follows that the third equation from (e 1 — If 
column must be a quadratic in 9 P ^{~ 1 , independent of & 1 , {= ^(^p* 1 "" 1 ) say}, and 
such that 






'j-1—0 n 

\j p w p 



Thus if p l x is one root of <j> 2 = 0, p l is the other. 
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Similarly, the fourth equation from the preceding column must give a cubic for 
P ~ 2 , (</> 3 = 0), such that 






dp 1 2 —0 p 

Thus 0/ 1 " 2 , P € ^\ p ei are the roots of <£ 3 = 0. 

Eventually the first column gives an equation of degree e x for P (viz., X ei ), of which 
the roots in any order are a possible set of values for P ...0 P \ Calling this equation 
(f) (#) = 0, and denoting its roots in some assigned order by cr u o- 2 , ...<r ei , let us consider 
the values determined for 1 p - l ... by taking P = ov.. and p l = <r ei . 

10. Again, as prior in order of simplicity, let the case in which the roots of (f)(0) 
are all different be taken first. 

It has been shown above that the equation Y„ when Xi... are substituted for y± . . . 
and 9 P , l p - 2 ... for #/, 2 p „ 2 ..., becomes identical with A^X^. 

Now, Y ei _! is merely a polynomial in P 2 9 independent of yi 1 ... and 2 p - x ..^ and 
vanishing for 2 = cr 2 , {r 8 ...cr ei . 

Let Y^ = ^(0/). 

Thus X ei+1 is linear in l p ~i the coefficient of the same being ^> ei _ 1 (^ 1 ); the part 
independent of l p ^ contains only #/, which is now a determinate quantity. 

If the roots of <f> tl (6 p ) are all unequal, ^ ei „ 1 (^ ? 1 )^=0 and ^_i is given uniquely; and 
similarly Y 6l = is a linear equation for 2 p - x in terms of #/ and #/, Z ei -x a linear 
equation in 0\- u and so on. 

It is important now to consider whether the order in which the roots cr 1 ...cr ei are 
taken is of significance in the solution ; that is, whether the value of k p ~ 1 associated 
with a particular root <r k is the same whichever column of the dependent variables 
this root is associated with, and whether a change in the order necessarily implies a 
distinct solution, because, if so, the solution would appear to be by no means unique. 

The equation X ei+1 giving X P ~ X is, we have seen, of the form <j) l (cri)0 l p -i-\-\Jj((Ti) = 0, 

where ?(ff) = ^(0). 

Thus 

ay - ^M-^^ ffl , MeAH&z) 

^A+i =: — u p-\\ * 

CTi — CT 2 CTi — CTa 

Now ^(o-x) is save for a constant factor 

(CTi - CT 2 ) (a*! - <r 3 ) . . . ((Tx - 0" ei ), 

so that A^X ei+1 is 

{(cr 1 -cr 3 )(o- 1 -<r 4 )...((r 1 -(T fi ) + (q 2 -cr 3 )...( ( T 2 -crJ}^- 1 + ^^^" ,j, ^ 2 > = 0. 

But the equation Y €l , which is independent of 2 p ~ 2 ,.„ becomes, when #/( = o"i) is 
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substituted for 6 P , the equation A^X^, so that when P is substituted for P in 
Y 6| it becomes symmetrical in 6 l p ^ and 6 2 p ~ x . Y €l must therefore be of the form 

A(<Ti, cr 2 ) ^_ 1 + B(cr 1 , cr 2 ) 0V-i + C(cri, cr 2 ) = 0, 

where 

A(cr 1? en) = B(cr 1? CTi) 

and 

A(o-i, or 3 ) + B((ri, cr 2 ) = (cTi — a 3 )...+(cr 2 — or 3 ).... 

XT 

Hence 

A((r b cr,) = (^-^(^-cr,)... = B(* u <r,). 

Now before the values of #/, P are substituted in Y ei the coefficient of 2 p - x must 
be a function of p only. 

Therefore B must be a function of o* 2 independent of ov 
Hence 

B = (ora-cr 3 )...(cr 2 -cr ei ) 
and 

A = (cri — ^...(o-! — cr ei ). 

The equation for fl 2 ^ is therefore 

a - ! — cr 2 
By virtue of the equation giving 0^-1, therefore, we have 

(cr 2 »cr 3 )...(c7 2 ~a- 6l )6Vi + ^ ( ^ = 

C7 2 <J*i 

as the equation for 2 p „ x . 

But this is identically the equation that would have been obtained for 1 p ^ 1 if 
a*!, cr 2 had been interchanged. Thus the value of P ^ X associated with the root cr 2 is 
unaltered by this interchange. 

We have further to see that the same permutation does not alter the P ^ X for the 
subsequent columns. 

In just the same way as above the equation for 0\^ x is shown to be 

(<r»-<r 4 ) (o- a -cr 5 )...^- 1 + (<7 8 ^ 
which gives 

(<r»- *<)... 6*^+ 7 ^ i = 

independently of the order of a- x and cr 2 . 

The same holds for k p - x > k < c l3 and a similar proof for the interchange of any other 
pair of the roots cr. 

Thus supposing the roots of <f>(9 p ) = to be all different, there is associated with 
each a unique determinate value of P „ X . 

VOL. CQV, — A, P 
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These quantities then being determined, consider now the equation X ej+2 , Y fi+1 ... 
and first it must be pointed out that the relations established between the equations 
X, Y... through the operators A (pp. 14-15), where the quantities pi p - ltt . were 
considered as independent, are still valid when p _ l9 &c, are determined as functions 
of P . 

The operator A p in the first place becomes replaced in the equation X e]+2 after that 

l _ — 6 2 _. 
giving l p - x by A^ + A^-!.--^- 1 - — Trip- ^ u ^ A^_iX ei+2 vanishes owing to the choice of 

u p u p 

l p -i> so that the value of 1 p ^. 1 being substituted in X e]+2 the operator A p may be still 

applied to establish a relation with Y €l+1 . 

We have further 

A p - 2 X ei+2 — (j) ei ^ l (0 p 1 ), 

while A^_ 3 X- i+2 ... vanish identically because of the vanishing of Y e] „ 2 , Y^.... 
Thus the equation X ei+2 is of the form 

in which p l and l p - x are to have their determined values, so that the equation may 
be written 

i r i(o-i)^-2+x(^i) = °' 

The operation A^ having been shown to be applicable to the equation in this form, 
reasoning exactly as above shows that the equation for &\- 2 reduces to 

so that the values of p - 2 associated with the roots cr l9 <t 2 are independent of the order 
in which these roots are taken, and likewise the values of # 3 ^_ 2 ... will be unaltered 
by a permutation of the same. The same may be shown in the same way of a 
permutation of any other two consecutive roots, viz., that such permutation gives rise 
to a corresponding permutation of the ^„ 2 .... Identical reasoning leads to an 
identical conclusion with regard to P ^ Q ...0 2 . 

Eventually we come to equations giving V When p 1 ...0 2 l have all been found, 
the equation X e]+j p_ r is of the form <j) ei - 1 (cr 1 )0 1 l + x( o "i) = 0> where, as before, the 
coefiicient of X is not zero, so that x l is determined like the rest; while 1 2 ...0 l * 1 are 
found respectively from Y ei+J0 _ 2 .... 

All the #'s being now determined, if we pass to the equation X ei+JP and follow the 
same argument that was required to prove the preceding equations independent of 
a?! 1 ..., the coefficient of x x is found to be the left-hand member of X ei+ ^_i with (#i x -f 1) 
substituted for x l > i.e., it is simply fc^^cri), which is not zero. Thus X, i+Jt? determines 
the first of the undetermined elements x^.... 

Similarly in X ei+J9+1 the coefficient of x 2 l is '2</> t:i - 1 (cr 1 ) 9 so that by this x 2 l is given, 
and so on for the succeeding equations in turn. 
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In order to proceed to the determination of the second column it may be noticed at 
once that the coefficient of y r x in any equation is exactly equal to that of x r l in the 
corresponding equation from the first column with o"i and a 2 interchanged, which 
includes the interchange of 6 l p ~ k and 6 2 p ~ h . In the equation Y tl+p - l9 therefore, the 
coefficient of y± is identically zero, while the unknown 6 2 s are now all determined. 
The closer consideration of this equation is deferred for a moment. 

The coefficient of y± in Y €i+P = coefficient of Xi in X €i+P with o^ and a 2 interchanged 

= (<r 3 — o-i) (0-2-0-3)... (<y 2 -<T € ) i= 0. 

Similarly the coefficient of y 2 l in Y €i+p+1 , and in general of y k l in Y €l+p+k -. u is not 
zero, while Y ei+p+k -i does not contain any element y v } for which m>&. 

For the third column the equation Z ei+j0 _ 3 determines 6* 9 and the two equations 
following, Z €l+p - 3 > Z tl + P -i> are still independent of z^ 9 z 2 ... 9 while the equation 
Z^+p+A-! contains z 1 1 ...z k 1 9 the coefficient of z k being &(cr 3 — cr 1 )(cr 3 — a- 2 )(a- 3 — cr 4 ).... 

In the same way, if the elements of the e^ 11 column be denoted by £ and the 
associated equations by ft, the equation Q, p determines 8i\ I2 1 , +1 ...flp +ei _ 1 do not contain 
& 1 ..., and i^ +ei contains & 1 only. 

11. So far the matrices x?---Xi ^ ave been taken to be simply diagonals. It will now 
be shown that the insertion of constants to the right of the diagonal in the first e l 
columns of Xi? can t> e carried out in such a manner as to affect none of the conclusions 
hitherto made, while they may be chosen so that the equations Y €i + p - ly Z ei+J9 _ 2 , 
Z^+p-x, &c, are all satisfied. 

Denoting by a {j a constant in the i th column and j th row, i>j> i^€ u j <ej, the 
Y equations become 

(0) (* P +i-0 P +i)yo-x o =o. 

(1) (ap+i-$ 2 >+i)yi+(a P -Qp)yo-Xi~oL 2l x ^o, 



(r+l) (^ P +i-O p+1 )y r+1 + (a p -0 p )y r ...-x r+1 -~a 2l x r =:O. 

These equations are to be treated just as they were before the constants a were 
introduced — the same elements remain undetermined as before, but at each stage the 
quantities found presumably contain a 21 . 

We see at once from equation (1) that the coefficient of a 2l in y x (the first element 
being excepted) is simply y . In fact, it can be shown step by step that the coefficient 
of a 21 in y r+1 (the first element always excepted) is exactly that part of y r which is 
independent of a 21 with 8 2 increased by 1 ; and therefore the coefficient of a 21 in Y r+1 , 
when the values of y f 2 ... as far as they are known are substituted, is equal to that 
part of the left-hand member of Y r which is independent of a 21 , 6 2 being increased by 
unity. 

Now Yx is independent of a 21 and of 6 2 and vanishes, thus Y 2 is independent of a 21 , 

D 2 
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and is therefore the same as if a 21 were zero. It is also independent of 2 , so that Y 3 
again is independent of a 21 . 

Thus until k is so large that Y* does not vanish independently of 6 2 -> Y k+1 is 
independent of a 21 , and therefore the same as was obtained in the foregoing, where a 21 
was neglected. 

Thus the insertion of a 21 in ^o does not affect any of the equations Y 1 ...Y ei+p _ 2j and 
therefore the values of 6 2 ...6 2 are independent of a 21 . 

But in the equation Y^+p^ the coefficient of a 21 is the left-hand member of Y ei+p - 2 
with X 2 +1 for 6 2 

= (<r 2 — cr 3 )(cr 2 — o- 4 )... ^ 0. 

Thus a 21 can certainly be chosen so that the equation Y^+p^ is satisfied. 

Having determined a 21 , it is at once seen from p. 19 that the following equations 
now determine y x l , y 2 --- without ambiguity; for since 6 2 ...6 2 are independent of a 21 , 
the coefficients of y*. &a, are those found there whether a 21 be zero or not. 

In the same way for the third column, with a 32 , a 31 , taken into account, the equations 
become 

\ a 2)+i~~~v p+i) %o y$ == "> 



and just as before the first equation in which a 32 occurs with a non-vanishing coefficient 
is the one following the equation from which X Z first does not vanish out identically, 
viz., Z, i+r _ 2 ; while a 31 will occur first in the equation homologous to the Y equation 
in which a 2l first occurred, viz., Z^^ ; in fact, in Z €i+i? _ 2 , a 32 will occur multiplied by 
the left-hand member of Z ei+J9 _ 3 with X Z +1 put for 0*, and in Z ei+j3 _ 2 , a 31 will be 
multiplied by the left-hand member of Y €i+P „ 2 with 0^+1 put for 2 : both these 
factors are other than zero. 

Thus a 32 can be chosen to satisfy Z 6]+;? _ 3 , and subsequently a 31 to satisfy Z €i+1 ,_ a , 
while the preceding equations are quite independent of them both ; just as for y, then, 
z x l ... are given in succession without ambiguity. 

Treating the remainder of the first e x columns in just the same way, all the elements 
of these columns are found in succession, and the solution is complete as far as these 
columns are concerned. 

The e 2 columns associated with the next group of equal roots may be treated in 
the same way, the singular equations being in this case the (ex + e^ 11 of each set; 
constants oC j will again be chosen in the matrix x ^° ^ e right of the diagonal, 
cx + cg^: t> Ci + 1, € 1 + € 3 >y > € x +l, to satisfy certain equations as above, and so for 
each root in succession. 

Thus if the various equations for P associated with the different groups of equal 
roots of the characteristic equation have their roots all different, and the " equations 
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of condition " (p. 15) for each root are satisfied, a formal solution of the linear system 
has been found in the form 



i«H?s+!f+-+¥). 



where the elements of yj are series of positive integral ascending powers of t, and 
Xi-.-Xp-i have all elements zero save those in the diagonals, which are made up of 
determinate constants ; and Xp consists merely of square matrices about its diagonal 
of € x , € 2 ... rows and columns respectively, each of which has zero everywhere to the 
left of the diagonal and determinate constants everywhere else. The elements of 77 
are in general divergent. 

The matrix O above will be known as the " determining matrizant." As occasion 
will be found later to discuss a more general matrizant, nothing further will be said 
of it here except for the case in which p = 1 , which will be worked out fully in order 
to make clear the march of ideas in the more general case. 

12. For p = 1 the equation X 6l is an equation of degree e x for X > Y €l ^ is of 
degree Ci — 1 for 2 > and so on. 

K X e| = * (<V) = 0, then Y ei . x = iM=M!=l} , 

so that the remaining roots are those of Y ei - X diminished by unity. 

Similarly the roots of Y ei -i are those of Z ei _ 2 diminished by unity, and so on, so that 
the roots of ^i 1 ) = are X \ 2 -l, ..., 0^-e^l. 

The equation 

X ei+1 is x 1 ^(0 l 1 +l) + x (0 l ') = O; 

x ei+2 is x^ (e 1 i +2)+x 1 \is (Ofi+x 1 (Q1 1 ) = o ? 

and so on. The equations for t/i..* are of the same form, with X 2 for ^ We 
suppose therefore in the first place that X , X 2 ... do not differ by positive integers or 
by zero, so that the coefficients of the first terms in these equations are all other than 
zero, and all the x s and y's are determined uniquely. The quantities a being then 
determined, as above, the solution is altogether determinate. 

If p = 1 and x l . . . 0{ l do not differ among themselves by integers, then the solution 
is of the form 



yjn< 



%\ 1 ...\ l /0 x l a 21 a a ... n 
OS 1 0...|p+ I X 2 * 32 ...J7 



v» 



_/ 



'& 1 



A 



'2 v "*" \ 1 / * "Ml 

ft 2 O.J^2 + 



^n< 



in which a ij = 0, unless #> l = &/ 



^ 2 O..J7 + 



* 




/0 a 21 ...\ 



a 



31 



\. 



1 
t 



r t tf t \ 



say, 
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Now £l(w + <r) = ft^ftjft^wjcrft^)} (Dr. Baker, Zoc. dt> p. 339), and 



e 2 i 



a 1^1+^) - o»ai -\ where <u = / e * ^ 



and w is the value of &> at £ = £ °, so that 







e»* 



c 



T/«i* 



V * * 2 * 



r 



= ^n. 



a; 



o 



1 

£ 2 






3 • * 



0, 






-0,2+0,3 



+ 



1 



0, a 



21 



t 0. 










a 



31 






<X 



32 / ^ \ 



\*. 



0' 



"\ 



v. 



e • 



• a 



* • o 



there being no exponentials in the last matrix since 

a ij = 0, unless #>* = 6J. 



e 2 l e\ (v\ 



Writing this v n ( § + El ) O ( 1 ), r is a matrix having zero in and to the left of the 



diagonal, so that Q ( — 2 ) has zeros in the same places. 

r /r\ 

~2 Qf-g) therefore has zeros in and to the left of the diagonal, and also in the (n—1) 
places to the right of the diagonal, and also wherever r has a zero, and so on. 

r r 

Thus Q70Q70... vanishes after a finite number of steps. Further, none of these 

t t 

expressions contain log t y since T contains no integral powers of t. 
Thus 

A, 721> 731 
'*p\ I 1 



... 1 



\ 



a 



t 



a 



733 

1 



* » 



* » 



all the places which were occupied by zeros in Y being also occupied by zeros in this, 
and jij contains only a finite number of powers of t, positive or negative, and no 
logarithms. 

We may specify a little more exactly the form of the term y {j . A typical term o{ 
Y/t 2 is 

t? \tj t \t , 



and 0,-Oi is not an integer and c {j is unity or zero. 
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The corresponding term in Q (T/t 2 ) is 



y 



A/ J L \&o-' 



r r 

It follows that in -§ Q -5 the r th column is a sum of terms belonging to the indices 

di—0i\ where s = 1, 2...r, and so for each operation Q. 
Thus finally we have the result — 

The term in the i th row and/ 11 column of &(~ 2 ) is a sum of terms belonging to the 

If) 0\ 
indices X S —0 X \ s = 1, 2...jf ; so that in the product O ( -f + — 1 the j th column is a sum 

of terms belonging to the indices 1 1 ...0 i j . 

13. Supposing still that p = 1, let the indices 1 1 ...0 x l cease to differ by other than 
integers, and likewise the other groups. Let them be arranged in groups differing 
by integers, so that their real parts are in descending order of magnitude in each 
group. 

Then no root of <£ (0^) = will exceed X by an integer, and therefore (j> (^i 1 4-^)^= — ^ 
a positive integer, so that the equations X 6]+1 , ... do not fail to determine x^.... 

If, however, 2 = l l —k, <f> (0 2j rk) = ; so that the coefficient of y k in Y ei+k vanishes, 
leaving y k unknown. We take y k = as the simplest assumption, and the following 
equations then give 2/Vn> &c., all without ambiguity. We are, however, left with 
Y ei and Y €i+k in general unsatisfied. Of these one can in general be satisfied without 
affecting the rest of the argument by an adjustment of the element x/ 1 . 

It has been seen that a constant a 21 in the matrix ^i occurs first with non- vanishing 
coefficient in Y £ . 

Clearly, then, if we introduce a?H ] \ it will leave all the equations to Y €i+k ^ unaffected, 
and add to Y €i+k the left-hand member of Y ei „x with di+k+1 for 2 . 

But Y ej _i is an equation of degree e l — 1 of which 2 is the greatest root, so that the 
multiple of a 21 added to Y €i+k is not zero. Thus a proper choice of a 21 satisfies Y ei+k . 

Again suppose X = 2 — k x = 0^—^ — ^ ^>0, & 2 >0. 

Then the equations Z €i+k9 Z €i+ki+kj fail to give z\ 9 z l ki+k /, but a 31 , a 21 can again be 
determined so that, if a sl t kl+k *, a S2 t kl occupy the places above 0* in Xi> the equations 
Z ei+Aj , Z €]+A . i+A . 2 are satisfied, 0* being unaffected and z x ki> z x ki+k9 being taken zero. 

Suppose then 1 1 ...0 1 h form the first group of l l ..,0 l €l differing by integers. Then 
treating the first h columns all in the same way, the ^h(h— l) equations Y €l , Z ei , 
Z €l _!... must be satisfied identically when the X have been all determined, and must 
be added to the equations of condition already found. 

Suppose now 1 h ~ {1 ...0 1 h+k form the next group of roots differing by integers and 
consider the (h-\-r) th column r<&. The equation giving x h+r is that indicated by the 
suffix €i— (h+r— 1), and those following this up to that with the suffix e x are 
independent of the undetermined elements of 7?. 
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Further 1 — 1 of the equations subsequent to these fail to determine the appropriate 
element as above, on account of the quantities 1 h ~ i ' s ~~0 l k+r being positive integers 
for s = 1, 2...(r— 1). These (r— 1) equations are satisfied by putting terms 

t e " ~ dl a h+r,h ' { ' k (k = 1...T—1) in )(i> while of the other h + 1 — 1 equations constants 
a A+r,s (s = 1...A) can be found to satisfy h. Thus (i — 1) equations of condition are 
found from this column, and therefore ^k(k-~ 1) from this group of roots; and so for 
each group of roots. 

Assuming all the equations of condition to be satisfied, we have now the following 
formal solution 



rjn 



a< 



Xi 



-J- £5 
* t 



where Xi m as follows : — 

The square matrices about the diagonal of h, /c... rows and columns respectively, 
corresponding to the groups of 0^... which differ by integers, are of the form 



0X 1 



21 • * * h 



e; 2 






i . 

i > 



and all other elements, to the right of the diagonal and within the matrices of e l ... 
rows and columns about the diagonal corresponding to the groups of equal roots for # 2 , 
are numerical constants, and all others to the left and right of the diagonal are zero. 
Applying now the formula 

n(w+a-) = a(w)a{a" 1 (w)orn(w)} 9 



the solution is put in the form 



r 



**(§ + 



Hi 

t 



D.A 



1 



/0, (t/t )-w, 



0. 



0. 



w 



■e^ + e? 



i + 



i 
t 



/o 



^21 * -*S1 • • • 



"^ 



Coo... 



J 



where in the last matrix all elements are zero that were zero in ^ 1; and Oy- is a 
constant if 0j—0i is a positive integer, but otherwise is a numerical multiple of t d) ~ e \ 

The expansion of the matrizant can be effected as on p. 22, with the result that in 
the expanded matrix the first h elements of the first row contain log (t/t ) to the 
powers 0, 1...A— 1 respectively, while the rest of the row is free from logarithms; the 
second row begins with zero, then unity, and the next (h—2) elements contain log (t/t ) 
to powers l...(h—2) respectively, and so on, the h th row being entirely free from 
logarithms. In the (h-\-l) th row in the k places beginning with the diagonal term 
log (t/t ) occurs to powers 0, !...(&— 1), and so on. 
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14. Returning now to the general case left out of consideration on p. 16, in which the 
roots of <f>(0 p ) = are not all unequal, we suppose the roots of this equation arranged 
in groups, of which the members are all equal. 

If n = cr 2 = ... = cr r9 ^-1(0-!) = fc-ifa) = 0, since the roots of <t> €l -i(0) = are 
o" 2 , ...cr r ; and again, <£ ei _ 2 (<r 2 ), ^ ei _ 3 (<7 8 )...^)( ei _ r+ i)(cr r _ 1 ) are all zero, where 
(f>^_ 2 = 0...<£ ei _ r+1 = are the equations for #/, P ..«0 P . 

The equation for l p „ x (p. 16) reduces therefore simply to xjjfa) = and fails to 
determine ^_! ; but, cr x being already known, this must be merely an " equation of 
condition " among the coefficients. 

Similarly in the equation X e]+2 the coefficient of l p „ 2 vanishes, and this equation 
therefore is of the form 

t/,(0V-i, #/) = or t//(^V-i^i) = 0. 

Now the operator A^_i acting on this equation, since cri = cr 29 gives the equation Y e , 
which, as has been seen (p. 17), is linear in 2 p -i 9 the coefficient being (£ ei _ 2 (cr 2 ). 

If r = 2, this does not vanish, and therefore X ei+2 is a quadratic for l p - X9 of which, 
owing to the relation through A^, 1 p ~ u 2 P -\ are the two roots. 

If, however, r>2, the equation Y ei must become an equation of condition, since 
<£ ei _ 2 (cr 2 ) = 0, and therefore also X ei+2 becomes independent of l p ~i and gives another 
equation of condition. 

Carrying on this reasoning step by step, we find that X ei+1 ...X ei+r - 1 are all 
independent of 1 p ~ 1 , 6 1 p -2--> while X ei+r is of degree r in 1 p ^ 1 and independent of 
1 p - 2 .... If any root of this equation be taken for 9 l p -i the equations Y ei ...Y e +r _ 3 are 
independent of 0%- l9 2 p ~ 2 ... 9 while Y ei+r _ 2 is an equation of degree r— 1, which, since 
it is derived from X ei+r by the operator A^_ 1? has for its roots the remaining roots 
01 j\. e ^+ r . 

Choosing one of these for 2 p ~ l9 Z e]+r _ 4 gives an equation of degree r— 2 whose roots 
are the remaining, and so on. 

Similarly, if <r r+l = ... = cr r+sy p - l r+1 ...0 p „ l r+s are given as the roots of an equation 
of degree s 9 and so for each group of equal roots or. 

Consider now one such group with the values of P - X r+1 9 ...0 P - X r+S obtained. 

Let the equations of which these are the roots be 

<M<w +i ) = o, h-iM-r*) = 0, ..., Me P -r s ) = 0. 

Then, if the roots of cj) s be all unequal, say = Ti...r„ 

<t> s ~i(ri) ^ 9 <f> s - 2 (T 2 )=fc0 9 ..., 
but 

<t>s-i(Tj) = 0, y^=i. 

The subsequent equations are then seen by the application of A 3 to be 
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^(r^e^r'+xsi^ oi) = o, 

<£ s - 2 (t 2 )^-/ + Hx s -i(t 2 , cr a , Ti, (Ti) = 0, 



r + l 



which, since the coefficients of p „ 2 do not vanish, at once give the values of p - 2 
these, as in the case of P - X when the <x's were unequal, can be shown, if the roots t x ... 
undergo a permutation, to undergo the same permutation, so that the same V ^ 2 is 
associated with any particular r in whatever place this r is taken. 

If, however, the roots r fall into groups of which the members are equal to one 
another, these equations again resolve themselves into equations of condition owing 
to the vanishing of ^> s _ 1 (r 1 ), &c. ; and, as before, the quantities p - 2 fall into corre- 
sponding groups given as the roots of equations of degrees equal to the numbers in 
the respective groups. 

The process may clearly be carried on as far as the determination of 6 2 by the use 
of the operators A^_ 2 . . . A 2 . 

A further remark should be made as to the finding of X , in connection with the 
operator A 1? which has been defined to be such that 

Suppose that 2 . . . $ 2 k are given by a set of equations 

where the affixes of the w's denote the degrees of the equations, and the roots of each 
equation are the remaining roots of the preceding after any one of them has been 
chosen. 

Suppose that of these 6 2 ...6 2 h are equal, so that 

Or 1 = r 2 ... = 0J\ r = p+l,p, ...2. 

Then 0,^(0/) = 0, *> k - 2 (6i) = 0, ..., ^^(0^) = 0, but <o k - h (d 2 k ) =£ 0. 
Then if the X equations following o) k be denoted successively by 

x o) k — 0, 2 co k = 0, ..., 

!&)/, is independent of Xi> &c, by virtue of o) k = and the preceding equations, and 
therefore 

so that ^ is also independent of ^ and must therefore vanish identically when 2 is 
determined. Hence also 2 o) k is independent of Xi 9 &c, and therefore 
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But in the same way 

A 1 ,(ift) A - 1 ) = (o k - 2 (6 2 2 ) = 0, 

so that i(o k -. 1 (0 2 2 ) must also vanish when # 2 2 is found and thus i<o k ^ 1 (0 2 1 ) = 0, so that 2 <a& 
is also independent of ^ x \ 
Ultimately we have 

<o k - h (0 2 h ) = A^o^+iW)}, 

= AiAx {2^-^+2(^2 )} = •••> 

and o) k ^ h {9 2 h ) =£= and is independent of 0^. 

Thus h,o) k {6i) is an equation of degree /& for 0/, since d 2 = 6J\ 0* = 0*.... 

Suppose its roots are v u v 2 , ... v h . 

Take 6± to be Vj. Then h -.iO) k -i is an equation of degree h— 1 for 0^, and its roots 
are (p + v 2 ) (p + v s )..* 9 as shown by effecting the operation A x . 

Choosing one of these again to be 9 2 9 say p + v 2 , h - 2 <*>k-2 is an equation for 1 8 9 whose 
roots are (2p + v 3 ) 5 (2p + i> 4 ).... Thus the quantities 01...0* are ^ 1? (p + ^ 2 ), (2p + ^ 3 )... 5 
(^-l)_p + %. 

In order to particularise the order of the roots v u v 2 , ..., they are arranged as soon 
as found as follows : — 

Let all those roots which differ from one another by integers be grouped con- 
secutively and let the arrangement in each group be such that r ~ 1 —0 r = or a 
positive integer. Suppose, now, the equation h <*>k{Qi) = is the equation X^. 

Then the coefficient oix-i in X^ +1 is ^(tf^+l) which, since no root of h o) k {6) = 
exceeds 6± by a positive integer, does not vanish. X g+i is moreover independent of 
x 29 ..., owing to the equations h -^k = 0, ..., being satisfied independently of #/. 

Thus X. q+1 determines x^. 

Similarly, X g+2 gives sc 2 \ the coefficient being ^(^-f 2), and so forth. 

Of the Y equations, that determining 6 2 is obtained from X^ by the operator A x . 
xb is, m iac u, jL q—<n» 

The coefficient of y± in Y q ~ p+1 is equal to the coefficient of x± in X^_ i?+1 with 6 2 9 
6 2 2 , ..., substituted for 1 1 9 Q 29 .... But X. q - P+1 vanishes identically as far as 6± is 
concerned and # 2 2 , ..., are the same as 9 2 , .... Thus Y q - P+1 is independent of y x \ 
Similarly, Y q „ p+2 ...Y q - 1 are all independent of the undetermined elements of y. 

Suppose now 61 —6 2 = X (a positive integer). 

Then h u k (6 2 +\) = 0. 

The equation Y q+k contains y± 9 ..., y k9 the coefficient of the last of these being 
htok(@i 2 +k) which vanishes for h = X. 

If, now, in the matrix ^ r (r = p, jp — 1...2) the second element of the first row be 
c r 21 , and in ^i be c 21 £\ the constants c will, as before ( p = 1), affect first the equations 
Y (? _ 1 , +1 ...Y (r . 1 and Y q+X respectively, entering into these with non-vanishing coefficients. 
Let c r 21 be determined then to satisfy the first p— I of these equations; Y q+1 ...Y q+K ^ 

E 2 
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then give ^...^V-i- Y g+A then fails to give y^> but c 21 can be chosen to satisfy the 
equation and y^ may be taken zero. 

The following equations then give the remaining elements in succession. 

This leaves the equation Y^ in general unsatisfied, and a further equation of 
condition is therefore necessary. 

Similarly, if 2 -~ X S = jx (an integer) of the equations Z, we can, by proper choice 
of c r 82 (r =^9... 2), satisfy the p — 1 following that which determines #i 3 , viz., Z q - 2p ; 
and just as a proper choice of the constants c 21 enabled us to satisfy Y q - p+lt ..Y q - u the 
constants c r dl can be chosen to satisfy Z g -. JP+1 ...Z g _ 1 . 

Thus two equations, Z g _^, Z ff , are left unsatisfied in general. The two remaining 
constants, c* 2 and c x 81 , are utilised to satisfy the equations Z g+/jt , Z g+A+/jt , in which the 
coefficients of zj~ and ^\ +ja vanish respectively. To do this the terms c 81 £ A+/A and &H^ P 
are inserted in the third column of Xi- 

If then the indices 01...0J be equal, or differ from one another by integers, exactly 
similar treatment applies for each of the first I columns of 77, the i th column furnishing 
(i~l) equations of condition. 

For the (Z+l) th column, however, 0/— ^ +1 , (r^h) is not equal to zero or a positive 
integer. Thus h <o k (0i l+l + m) does not vanish for any value of m, and the Ip constants 
c r l * lt8 (r =jp...l, s = I... I) can be determined to satisfy the Ip equations between 
TJ g -ip and U^ +1 , XJ denoting an equation of the l + l th column, and, in particular, \J q -i p 
being the equation determining 0i +l and U g+1 determining u x l . 

For the next column, however, l l+l —0 1 1 * 2 may be a positive integer, X 1 say, so 
that h o) k (0 1 h+2 + X 1 ) = 0, and the (q + A 1 )* 11 equation, instead of determining the 
appropriate element of 77, can only be satisfied at the expense of the q tYl , by making 
the element above 0i +2 in j^— c l+2,l+1 t x \ The q th equation then becomes a further 
equation of condition. Thus we shall obtain r—1 equations of condition from the 
(l + r) th column, associated with an index belonging to the second group of indices X 
differing among themselves by positive integers ; and so on through all the indices as 
iar as U\ . 

A similar treatment is now applied to the columns (A +1 )...(/?, + &), for which 
r h+1 = #/ i+2 ... = 0J l+k (r = jp...2) ; x h+l is given as the root of an equation of degree h ; 
and the minors of the determinants ^; r , whose diagonals are r h+l ...0 r h+k , have the 
elements to the right of the diagonal suitably adjusted as above, while one equation 
of condition is furnished in connection with every difference 1 h+r —0 1 h+ % which is an 
integer. 

Supposing these equations all satisfied, consider the expansion of the matrix 

o(^~+... + ~), which is effected in just the same way as for p = 1 (p. 22). 

If w = 2 -£, where r is a matrix made up simply of the diagonal terms r l ...0 r n , 
the application of the equation 
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XI {w + <r) = fl (w) O [XI * (w) crXl (w)] 

resolves the required matrix into a product of two, of which the first, O (w), consists 
simply of diagonal terms of which the r th is 



r 



e 



_ d p+i 



6 



r 



3L 



0, 



r 



(p~l)t p 



-1 



t r 
t 1 



e' 



p±1 



e 



W 



f i 

^0 



The second matrix has zero everywhere in and to the left of the diagonal ; and, 
since in the matrices x "the element in the r th row and s ih column was zero unless 
6 J = m s (m = p+1, ...2), it contains no exponential expressions. It can therefore be 
completely integrated in finite terms, just as was done in the case p = 1 (p. 23). 

15. A simple example may be appended of the application of the method to a particular system. 
Consider the equation 

xy + xhj f (2 -x) +3/(1 - 3x + x 2 ) = 0. 

Putting y x ^y,y 2 = xhj\ we have the linear system 



¥ = 




y 



i_ 

X 2 

fcV d/J 



2/a; 2 V3 "' 







3/ a; 



1 



N 
0, 



y 



The characteristic equation is 



giving equal roots - 1 for p. 
1 CT 






1 
2-p 



-0 or ( P + 1)2 = 0, 



With p 



1 1 



the transformed system becomes 



y 



= uy say. 
Considering now the subsidiary equations 



/ - 1 1 \ 1 t /0 0\ 1 t / 

+ \0 3/^-1 



1/a; 2 



s 




y 



?/ ~U7] — rj 

Then the equations to be satisfied are 



-1 1 U . M Ml 
~W 10 2 /a? 



2/ # 



, and let w = ( x ^ l 

& y% 



I. ^i-flx^O. 
II. (1) a 2 2 - (1 + *i) ^i 1 = % (2 - 0i) ^ - 1 

I. and IL (2) give (fli- 1) 3 = 0. 

We take 6\ - 1 therefore, so that xi = 1 and £ 2 2 = 2a?! 1 . 

Again 

III. (1) £ 2 * - 6>! a**- 0, (2) (1 - X ) x£ - x? = 0, 

of which (2) gives Xi 1 = 0, so that x 2 2 = 0» 



0, 
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Similarly Xi 2 — x 2 B = 0, and so on. 

The equations for the second column are 

I. (1) gfe 1 - ai 1 - * = 0, (2) 3-0 2 -ai 2 = O, 

(2) gives 0.2 = 2 — 1 + 01 and (1) gives y^ 1 = A. 

II. (1) y£ - (1 + 2 ) yi i - ^ - X a** = 0, (2) (2 - 9 2 ) y£ - x? - \x£ = 0, 
(2) gives X — and (1) gives y 2 2 = 22/1 1 and also ^i 1 = 0. 

III. (1) ^ - % X 2 - ^ - Xx, 2 = 0, (2) (1 - 2 ) ^ - yi l = 0, 

so that yi l = and ^ 2 2 = 0. 

It is easily seen that all the remaining terms vanish. 
Thus the solution reducing at x =■- $o to the matrix unity is 

y \x l]\x*\Q -l/ + iV0 2/jU 1, 



where 



-G ?H<^(2 SM- 1 . * 



u=oi" 



L/-i o\ 1/1 o 
\x*\ o -1/ x\o 2/ r i 




6 x oXq \ Xq* 



Thus 



i 







e^oafo \ av I. \ a; 

1 Q\e^x I 1 \ o /0 — \/ l 0' 
„ 1 I" (o l)°( H( ,o 1. 





i 



1 0\ A 1 ! x 



1 (l IlogS' 

e-*a? /I 0\| £o ^o 

1 U 1/L 



a? 




16. The number of conditions found in the course of the analysis shows that the 
solution in this form — which may be called the " normal" form, by analogy with the 
name ' ' normal integrals " of linear equations — is by no means always possible. As 
many writers have pointed out, there is a much more general type of solution than 
the normal series for the ordinary linear equation, in the form of a normal series in a 
new independent variable x vk , h being a positive integer (Cayley, Hamburger, 
Fabry, &c). 

The method developed in the foregoing is peculiarly adaptable to the investigation 
of these integrals, inasmuch as the transformation to a new independent variable is 
very simply effected. 
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If in the linear system 

-7 = u(t)y we put t = <£(#), 

we have without any calculation the new system for y as a function of z 

Jl=<j> f (z){u(<$)}y. 
Suppose now <f>(z) = z k . Then the transformed system is 



-JL = kz k x u (z k ) y. 



y 



dz 
If, then, the original system is of rank p, so that 

the new system is 

dy hzp+x k(t v -uJra«z k - l + 

and is of rank &p. 

If, now, we were to put z = & 1//cp . 2 1 , the form of the equation would be unchanged, 
while the coefficient of z llkp+1 in the right-hand member would become the original 
canonical matrix a p+1 . This is, however, not necessary, as the whole investigation 
could be carried out equally well if any constants whatever replaced the unities to 
the right of the diagonal in ol p+1 . 

It may now well happen that though all the equations of condition found for the 
general system are not satisfied, those associated with the new system are all satisfied, 
so that the new system possesses a solution in normal form. If this is so, the original 
system will be said to admit of a solution in subnormal form. In fact, an integer k 
can always be found such that this is so, owing to the vanishing of the coefficients of 
z~ kp+r {r= 0, l...(£-2)}. 

In the first place, all the conditions arrived at from the equations X 1? Y l9 ... will 
be satisfied (p. 13), for the coefficient of z~ kv is identically zero; in general, the 
left-hand members of X n Y r . . . are rational integral functions of the elements of the 
matrices A^, A^„ l5 ..., A^_ r+1 , if A m stands for the matrix multiplying z~ n \ and 
contain no term independent of these elements. 

Now the conditions found on p. 13 arise from the equations X^.X^-!, Y x ...Y ei _ 2 , ..., 
and therefore involve the matrices A kp , ..., A^„ ei+2 . These conditions will therefore 
all be satisfied if Jc >: e^ Similarly, the analogous conditions for the second group of 
equal terms in the diagonal of a p+1 will be satisfied if k >: e 2 , and so on. 

Consider first, as being simplest of explanation and as containing the essential 
features, the case in which all the roots of the characteristic equation are equal, so 
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that €j = n. It will be shown that a subnormal form satisfying the equation certainly 
exists if a p ln =/= with h = n. 

We know from the foregoing investigation that 6 l np is given as the root of an 
equation of degree n, and that, if the roots of this equation are all different, no more 
conditions than those just mentioned as satisfied are necessary to ensure the existence 
of the subnormal solution. But in this case the equation for 6\ p is particularly easy 
to construct. We have, in fact, 

ll/JU^ "~ ~~ V nm — VJ, t/-/J OVJ — • • • — — t<(/J — VJ , 

IuJLq ~~ ~~ t' jin^X ~~~~ ^ J t *- / 2 ~~* * * * — 2 "" — ^ J 

/M 'Y* ^ . /?* 'Y* ^ . ffi- <W ^ - — /y» 5 . — — /y» ft , — A 



/f/^o w _i v n p*Aj n —2 u np ,v n—S ••• — u j 

_#1 ™ M -4- \ ln — 

from which at once we have 

1 

\Q up) , Al» _ a Ar m _ '^\A w n(p _ 1) + 1 } w 

Thus, unless A ln n{p - l)+x = 0, the values of # x ? ^ are all different, and a solution in 
subnormal form is therefore possible, as stated above, with the independent variable 
changed to x lln . 

If, however, A 1 \ ( ^_ 1)+1 = 0, we have 

1 = = ft n = 

and it will be found that the same conditions are necessary between the constants 
An(jp — 1) + 1 as were found previously (p. 13) between the constants a p9 A n(i? _ 1)+1 
being the same as n. a p , e.g., from the equations 

we have 

A 2, \ ( ^_i )+ i + A 1 ' w ~ 1 w(2? - 1)+1 = 0, ^'.e, d p 2,n + oL p ltn ~' 1 = 0, and so on. 

Consider now what happens when these conditions are not all satisfied. Suppose, 

for instance, cL p 2,n + ap' n ~ 1 =jfc= 0. Let the original system be transformed by the change 

of variables 

x~z k , h = n—1. 
Then from the equations 

h . z x — u kp = 0, 

fc . Z2 f j^i :=: u , ... 

Z. ~n /31 ^ ra-1 1 AU-1 a /31 „n i A 2 i» /y 2 C\ 

& • % n-l u kpZn-2 "T-ci- 7c(^-l) + l — u ? ^ ^ w-1 » ■**■ &(p-l)+l^l — ^> 

we obtain the equation for Q x kp 



SATISFYING LINEAR DIFFERENTIAL EQUATIONS. 33 



(Vkp) n _v kp / Ai, n -1 -4-A 2 ' n ^ — 



& w & 



) 



and since the last written coefficient is not zero, the roots of this equation are all 
different, and therefore the transformed system possesses normal integrals free from 
logarithms. 

Again, suppose a p hn = and that a 1 , 1,B "" 1 + a i , att also vanishes; then the roots of this 
equation all become zero, and we find that the condition a p ljn ~ 2 + a p 2 ' n ~ l + a p z,n = is 
also necessary. If this is not satisfied and the original system be transformed by 
x = z k 9 k = n— 2, the equation for 6\ p becomes 

i — [ -J— {J\. p ( k -. 2 ) + l-r J *- p(k~2) + l'-T J ^- p(A-2) + l/ ~ U ? 

which has two zero roots and the rest all different. If the one condition necessitated 
by the equality of the two zero roots is satisfied, the solution is again found. If this 
condition is not satisfied, the transformation z = £ 2 effects what is required. 

Suppose now all the conditions of p. 13 are satisfied. Then whatever value of k be 
taken, we have 6 r kp = 0, so that, as far as we have seen, the transformation does not 
render the solution any nearer. 

We must, in fact, proceed to consider the further conditions for the case when the 
roots''^ are equal (p. 13). 

Suppose, for instance, the first of these conditions is not satisfied, then putting 
k = ft, we shall have 6 l kp = . . . = 6 n kp , the conditions then necessary before the deter- 
mination of kp ~i will be satisfied, and we shall eventually obtain a binomial equation 
for Q kp - X of degree n in which the constant term does not vanish ; the roots of this- 
equation being all different, the subnormal integral exists. 

Thus we may go through all the equations of condition in turn. 

In the more general case, where the roots of the equation for 6 P fall into more than 
one group of equal roots, the procedure is exactly similar. 

For example, suppose that a/' 61 , a/ ,+1,ei+e »; ... are all different from zero. Then the 
preceding work suggests that a solution may be found in which the first e x rowa 
proceed according to powers of t Ve \ the next e 2 according to powers oft 1/e * 9 and so on. 

The whole would thus be of normal form with the variable t ljk 9 k being the least 
common multiple of ci, e 2 .... 

In fact, if we change the independent variable to t 1/k , k having this value, the matrix 
Aty_ €l+ i is identically zero, and the indices #V-#% are the roots of — (0 kp Y l + = 0. 

These roots are all equal, and the corresponding equations of condition are all 
satisfied owing to the vanishing of the matrices other than A kr + L . 

If, now, we form the equation for 6 l kp -. x we have 

/cx 2 " kp ~\ = 0, KX± u kp ~\X 2 = 0, . ...... -»•"♦, 



kp ~iX 2ex - 2 + A €l kp ~ 2si +i — 0, 
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giving 



f)l _ 7, ; A lei 



£j9 — 2e 1 + l/& 



£ 



If — = 2, A lw A y_ 2ei +i = k . a p , and therefore if o^ lei =£= 0, the € x roots of this equation 

are all different. If, however, - > 2, A^_ 2 e 1+ i again is identically zero, and the 

necessary equations of condition are again satisfied. 

Proceeding thus we find, in fact, that if a p Ul y£ 0, 6 s kp , ..., $ kp - k!ei+2 l s = l...€i all 
vanish and that 6 s kp -. k/ei+l is the root of a binomial equation of degree e 1? whose roots 
are all different, and so for the other divisors of k If, however, a p Ul = we have the 
same equations of condition again necessary, viz., a p 2ei and a p l > ei ~ 1 = 0, &c. 

Assuming, then, that a p Ul zfc 0, we find, without difficulty, that all the quantities 
r * vanish for s = l...e x , save those for which r is of the form —k^p—mje^ + l) so that 
the exponential arising in the first e x rows involves only t 11 " 1 and not t ljk . 

The discussion of whether the solution of the subsidiary equation proceeds according 
to powers of t l/ei only in the first e x columns will not be carried out in full here. It is 
enough to know that, provided a/' ei , a p 1+1,ei+e2 do not vanish, a subnormal form 
certainly exists satisfying the equation. 

If, however, one or more of these quantities does vanish, and one of the consequent 
equations of condition is not satisfied, we may, as on pp. 31-33, find a new integer 1% 
such that the necessary conditions for the existence of the subnormal form are 
satisfied. 

17. As a concluding example consider the system derived from the equation of third order and rank one 

,/// , %o£ 2 + a 2l z + a 22 , a ^ + a 3 iz 2 + ^32^ + 033 , . 

y + y + . -_ y-Kj, 



which with 
gives 



yi =y, y-i= »Y, ys = «y ', 



r 






^33 



1 
1 

a 2 2 0/ 






0\ 





2 

— <^32 ~ (X 2 \ 4 , 



-8 









-^30 




The characteristic equation is - p d - pa 2 2 + ^33 = 0. 

We shall confine ourselves to the case in which this equation has three equal roots. These must all 
then be zero, and a 22 — 0, a s $ — 0. 

For the equation then to possess a normal solution we must have a 2 \ = 0, a 32 = 0. Supposing these 
conditions not satisfied, put z = t 3 ; then the equation becomes 



y = 3 



i 



1 0\ 

001 -+-+-+ 
i* P P ■ 

\0 0/ 



\ 






%2 





2 







^ 



1 



i|+... >y. 



The subsidiary equations then become 
3^2 _ 3 l - 0, x^ = • 



3% 2 - 3 V - ftfo 1 ^ <V = 0, 3^ 3 3 - O^xJ - 2 W 



0, 



W - ^W - 3a 32 - 0. 
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The last gives - -H^ 1 ) 8 " 3a 32 = 0, so that ^ 3 1 = - 3 (a 32 )*, where any cube root may be taken, the other 
roots giving <9 3 2 , 3 3 . 
Further 

3% 2 - <w - <W - (i + ^i 1 ) ^i 1 = o, 

3% 3 - 6 % W _ ^1^2 + (J - 0! 1 ) % 2 = 0, 

- #sW - ^2 3 + (3 - 01 1 ) xj - S^^i 1 - 3o2iaJi 8 = ; 

so that ^ 2 1 = (fail (am)* > a nd # 2 2 , # 2 3 are given by taking the other roots for (a 32 )i 

Lastly the equation 

- # 3 W - 2 W + (2 - <V) x£ - 3fl 8 2«2 1 - 3« 2 i« 2 2 = 



of which the last gives 



gives 



0,i /^ % 3 + ^! %2 + \^t ^ - 0,1 0£ c^ + ^ ^\ + (2 - V) fe* - 308^1 - 3a 21 ^ = 0, 



which gives 

so that ^i 1 = 5. Similarly, #i 2 = #i 3 — 5, and a subnormal form, exists satisfying the equation, of which the 

2 .1 A. 

first column has the determining factor e 2z * a *** z * £f, and the other columns have the same factor with the 
other cube roots of a 32 . 

We may remark that this agrees with the results obtained for this equation by the ordinary methods 
(Forsyth, "Linear Differential Equations," § 99) under the assumption that a 32 ^0. We have shown 
that this is a necessary condition for the existence of the subnormal form in the variable t = zh satisfying 
the equation formally, unless we have also a 2 i = 0. 

If, however, « 32 = and a 2 i^0, then, as we have seen above, the transformation t = z§ will give us a 
system admitting of 3 normal solutions; the equation for 9\ is, in fact, (0{f - 4ft 2 i#i = 0, giving 6 1 = Q or 

± 2a 2 i*. 
We see, in fact, that, when a 32 = and a 2 i i* 0, the characteristic index of the original equation is 2, so 

that there will be one regular form satisfying the equation, i.e., an expression of the form af P (x). 

If a 32 , C&21 are both zero, the equation is of Fuchsian type. Thus the normal or subnormal forms are 
found in all cases. 



F 2 



